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Overview 



• SOFIA Overview 

• The Thermostructural Concern 

• Determination of Governing Parameters 

• FEM Model Development 

• Results 

• Conclusions 
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Stratospheric Observatory For Infrared Astronom 



• Highly Modified Boeing 747-SP 


• 17-ton Infrared Telescope: 




Telescope Cavity 
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- Primary Mirror: 2.5m diameter 

- Optimized for infrared 
wavelengths that cannot be 
accessed by any ground 
telescope or current space 
telescope 

Max Opening (shown): 58° 
Mobile observatory platform 


(anywhere, anytime) 

expansion complete-. : • 
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science flights begun 
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The Thermostructural Concern 


• Primary: Original SOFIA design included telescope cavity ground pre-cool, 
and for various reasons needed to consider flight testing without 

- Would door opening at altitude result in a thermal-shock strong enough to 
damage imager optics? 

• Secondary: Some parties concerned that flight test with original design 
would still have unsafe thermostructural loading due to air temperature 
change along flight path (will not fly isothermal flight path) 

• Tertiary: CTE mismatch, already mitigated by ground testing of imagers 

• 4 different optical components identified 

- FFI (Fine Field Imager) 

• Schmidt Plate (higher concern) 

• Achromat 

- WFI (Wide Field Imager) 

• Achromat 1 (higher concern) 

• Achromat 2 

• FFI Schmidt Plate analyzed, none others due to results of analyses and 
different fixture comparison (clamps vs. fixing rings) 
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Telescope Assembly 
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FFI Front Group 



TFAWS 2011 - August 15-19, 2011 


Determination of Governing Parameters 

• Determination of physics to be modeled 

- k Convection-induced thermal gradient (thermal stress) 

- ik Bolt pre-load induced stress (and CTE mismatch +/-) 

- ★ Circumferential clamping (CTE mismatch) 

- O Vibratory stress 

- 9 Acoustic pressure 

- O Max flight load 

• Determination of domain 

- Relatively thermally isolated front group containing Schmidt 
Plate was clear choice for geometry 




*Already mitigated by ground test of imager 
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Determination of Governing Parameters 



• Determination of model key result(s) & acceptance criteria 

- Glass is much stronger in compression than tension 

- Tensile strength strongly dependent upon surface finish 

- Glass manufacturer (Schott) TIE-33 “Design strength of optical glass 


and ZERODUR” 


• Low stress level allowable for infinite part life (8 MPa and 10 MPa allowable 
for optical & ZERODUR, respectively) 

• When application requires higher stresses, statistical approach provided, 
characteristic strength values for zero failure probability all > 20 MPa even 
for the worst surface finishes - this part is optically finished 

• So o p1 < 8-1 0 MPa adopted as target value rather than requirement based 
upon strength data available, but 20 MPa would be seen as conservative 
benchmark value for margin determination 

• These are low values: model aggressively conservative and if no positive 
margin then refine conservativism 

- Due to clamping mechanism & ground test, in addition to determining stress level 
it was important to determine stress composition 

- Because of uncertainty in some inputs, a parameter sensitivity study would have 
to be performed to make this analysis complete 
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Determination of Governing Parameters 



• Determination of material properties 

- Materials known, no batch properties, combination of 
manufacturer supplied and database, if “pencil sharpening” 
required then examine more closely (along with other layers of 
conservativism) 


• Determination of structural loading 

- Loading due to bolt pre-load in clamping mechanism 

- Agreed upon value of bolt pre-load range 2-3 kN (3 kN value 
seems strongly conservative) 


• Determination of thermal loading 

- Instead of working up cooling rate range, assume limiting 
scenario of convective thermal shock 


- 3 Governing parameters for convective thermal shock: initial 
temperature (T,), fluid temperature (T R ), convection coefficients 
(h) 


TFAWS 2011 - August 15-19, 2011 


9 



Determination of Governing Parameters 



• Determination of thermal loading (cont’d) 

- Initial temperature: 20 °C reasonable assumption given onboard aircraft 
systems 

- Fluid temperature: 

• Airflow ingested into telescope cavity, into FFI Baffle, impinging/over 
external surface of Schmidt Plate 


• Not freestream temperature, but recovery temperature of fluid 

- Determined using the isentropic, subsonic compressible flow equation, but 
modified to assume non-zero flow at Schmidt Plate surface 



i+r^K-mJ 2 



Where T R is the recovery temperature, R is a factor (0.9) to compensate for the 
process not being perfectly adiabatic, IVL and "L are freestream values, and M res 
is the residual flow velocity (a max value for the whole cavity being =0.1 ) 

• This leads to a conservative value of T R = -40 °C for max door-opening altitude 

• The resulting shock value (T, - T R ) = 60 °C 

• This is conservatively the worst possible scenario, there can only be less severe than 
this (finite air temperature cooling rate, smaller AT, etc.) 
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Determination of Governing Parameters 



• Determination of thermal loading (cont’d) 

- Convection coefficient 


• Dependent upon geometry and flowfield 

- Used CFD results for velocity range 

- Calculated several correlations, subsonic stagnation point @ 15 kft 
(lowest door-opening altitude) was the highest, used conservative 
velocity, h = 60 W/m 2 K 

- Flow around the body behind the headring (low speed/free, h = 5 
W/m 2 K) 

- It should be noted that a physically impossible AT & h 
combination (from different altitudes) leads to a very 
conservative analysis 
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FEM Model Development 

• Began with hand calc a = E a AT (35 MPa) for bounding value for model results 

• 2D structural model to look at mesh density for solution convergence to this 
idealization 

• 2D & 3D transient thermal models to investigate thermal response between 
components & t max grad 

• 3D coupled thermal-structural transient FEM with structural and thermal contact (to 
allow DOFs between plate and fixture) 

- Mesh refinement to allow contact calculations to run without physically impossible load 
concentrations 

• Working model -16,000 hex elements using 21,800 nodes, -1 day runtime 

• Model used to iterate on 4 key input parameters to investigate sensitivity to 
uncertainty 

- Clamp pre-load 

- Friction coefficients 

- Convection coefficients 

- Shock strength 

• Epiphany: another mesh refinement (in contact region) to determine convergence 
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FEM Model Development 


Coarse Working Mesh 
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FEM Model Development 

MSC.Patran pre/post processor, MSC.Marc solver 
Loading separated for stress composition determination 



Non-linear static 

Initial State 

Non-linear transient 
load step (clamp pre- 

ioaa step (ciamp 
pre-load only) 


load + thermal 
loading) 
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Results: Max Principal Variation with Pre-Load 
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Very conservative (physically impossible) combination of AT = 60 °C, and h 
= 60 W/m 2 K produced only small increase in maximum occurring o p1 over 
clamp pre-load induced level (pre-load > 90% of max occurring) 
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Results: Max Principal Variation with Friction 



Time (sec) 



• 3 kN pre-load, h = 60 W/m 2 K, AT = 60 °C 

• Friction coefficients variation did not effect Schmidt Plate stress state 

• High-p case also run (not plotted) and produced overlapping results 
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Results: Max Principal Variation with Shock Strength 



Time (sec) 


• 2 kN pre-load, h = 60 W/m 2 K, AT = 30/60/90 °C 

• For AT = 0-90 °C o p1 only 1 .5 MPa higher 

• This would be less given a more realistic convection coefficient 
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Results: Max Principal Variation with Convection Coefficien 



80 W/m2K 

60 W/m2K 

40 W/m2K 


Time (sec) 

2 kN pre-load, AT = 90 °C shock, h = 40/60/80 W/m 2 K 

Higher shock value used to magnify the insensitivity 

For h = 0-80 W/m 2 K o p1 only 1 .5 MPa higher with overly conservative 90 °C 
shock 
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Results: Additional Refinement 


• Until this point several mesh refinements had been 
performed (2D thermal, 2D thermostructural, 3D thermal) 

• Additionally, there was a mesh densification required in 
contact regions to protect against false stress 
concentrations 

• The preceding two steps engendered a sense of 
sufficiency in model development, convergence 

• Upon realization that solution convergence needed to be 
demonstrated with “working” mesh another mesh 
refinement was done 

• Edge length reduction of 40% lead to max o p1 dropping 
from 10.6 MPa to 0.7 MPa, with further edge length 
reduction providing negligible reduction, proportional 
reduction in variation due to thermal effects 
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Conclusions 

• Determined clamp bolt pre-load contributed >90% of 
max principal stress 

• Demonstrated thermally-induced stress variation as 
insensitive to uncertainty 

• Found positive margin for extreme scenarios, no fatigue 
concerns 

• Cleared, by comparison, other components 

• Learned valuable lesson - keep careful track of steps 
taken (and any remaining) to demonstrate solution 
convergence when performing multi-disciplinary 
analyses, regardless of whether using home-grown or 
commercial code 

- Sense of conventional mesh sufficiency may no longer apply 
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Questions? 
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